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There are two complementary approaches to reaUzing quantum information so that it is 
protected from a given set of error operators. Both involve encoding information by means 
^ - of subsystems. One is initialization-based error protection, which involves a quantum op- 

■ eration that is applied before error events occur The other is operator quantum error 

O . correction, which uses a recovery operation applied after the errors ||2|]. Together, the two 

^ ! approaches make it clear how quantum information can be stored at all stages of a process 

I involving alternating error and quantum operations. In particular, there is always a sub- 

system that faithfully represents the desired quantum information. We give a definition of 
^ ■ faithful realization of quantum information and show that it always involves subsystems. 

CN . This justifies the "subsystems principle" for realizing quantum information fs]]. In the pres- 

ence of errors, one can make use of noiseless, (initialization) protectable, or error-correcting 
J> I subsystems. We give an exphcit algorithm for finding optimal noiseless subsystems by re- 

^ ' fining the strategy given in 0]. Finding optimal protectable or error-correcting subsystems 

^ ■ is in general difficult. Verifying that a subsystem is error-correcting involves only linear 

2 ■ algebra We discuss the verification problem for protectable subsystems and reduce 

\^ , it to a simpler version of the problem of finding enw-detecting codes. 

o 

a,: 



PACS numbers: 03.67.Pp, 03.67.Hk, 03.67.Lx 



^3 : I. INTRODUCTION 



According to quantum information theory, quantum information is represented by states of a 
^ ■ number of qudits, which are idealized d-level quantum systems. Here, we consider finite quantum 
information and do not explicitly refer to the underlying tensor-product structure of the state space 
of more than one qudit. Thus, we consider quantum information represented by the states of an 
ideal system I whose state space is determined by a Hilbert space 'H\ of finite dimension N. In 
order to realize quantum information in a physical system P, two problems need to be solved. 
The first is to determine the ways in which the state space of P can usefully encode the desired 
quantum information, the second is to determine which among these ways can be best protected 
from decoherence due to the dynamics of P and its interactions. If quantum information is intended 
for use in quantum algorithms, a third problem is to ensure that the dynamics of the system and its 
interactions can be used to implement the desired quantum gates. Here we consider the first two 
problems. The mathematical conventions for objects and notations are explained at the end of the 
Introduction. 
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A. Encoding Quantum Information and the Subsystems Principle 

A general method for realizing or encoding quantum information is as a subsystem of Tip Jsl]. 
This method involves a decomposition of as 

Up = {Hv ® Hs) © Hr, (1) 

where © denotes the orthogonal sum of Hilbert spaces. This decomposition identifies I' as a sub- 
system of P, written as I' ^ P. We use "primes" (as in I') for systems with identical state spaces. 
With this convention, Eq. ^ yields a representation of the states of the ideal quantum information 
system I in a subsystem of Hp. We say that I is "encoded" in P and call the decomposition of 
Eq. dl]) a "subsystem encoding". For lack of a better word, we refer to S as the "cosubsystem" 
of I' in P. R is the "remainder system". An example of such an encoding is that of a vibrational 
qubit of a single ion trapped in a one-dimensional harmonic potential. In this case the state space 
is spanned by the internal and vibrational levels of the ion. The space 7i|/ © Hs is formed from the 
first two vibrational levels and the internal levels of the ion, respectively. The space Hr consists 
of states with more than two vibrational quanta irrespective of the internal state of the ion. The 
familiar cases of such encodings have the property that H\i belongs to a physically meaningful 
degree of freedom. However, for the purpose of protecting against errors, it is usually necessary to 
use "entangled" encodings. An example is the noiseless qubit encoded in three spin-1/2 particles 
subject to collective decoherence [1]. A feature of subsystem encoding is that the states of I are 
not uniquely encoded as states of P. This is because any change of state of the cosubsystem S does 
not affect the states of I'. 

Subsystem encodings of I in P are equivalent to f-preserving isomorphic embeddings of B{H\) 
into 5(7^p) iHH. Here BiTL) denotes the set of (bounded) operators on 7-^. In particular, given 
the subsystem encoding of Eq. ([l]), B(H\) is isomorphic to the algebra of operators of the form 
A © I © with A acting on Hy, I on Hs and on Hr, where the operators are transported to 
Hp via the isomorphism implicit in the subsystem decomposition as needed. Conversely, if A is 
a subalgebra of bounded operators on Hp and A is f-isomorphic to B(H\), then there is a unique 
subsystem encoding such that the operators of A are the operators of the form A © I © as above. 

Are there ways of encoding quantum information that do not involve a subsystem encoding? In 
an attempt to answer this question, it is worth considering other prescriptions for encoding quan- 
tum information. There are two operationally defined ways of characterizing encoded information. 
The first is by a traditional encoding operation that isometrically embeds H\ into Hp. (Isometrics 
are linear maps preserving the inner product.) This is the prescription used in the traditional theory 
of quantum error correction and corresponds to a subsystem encoding with trivial cosubsystem. 
In this case, the subspace Hy © Hs = Hy of Hp is known as a "quantum code". That this is 
inadequate is apparent when one considers enlarging P by other relevant degrees of freedom. Sub- 
spaces also fail to capture the location of quantum information in realistic error-control settings, 
in particular fault-tolerant quantum computing. This is because in practice, error control never 
results in restoration of the encoded quantum information to any fixed quantum code. Assuming 
that this is a requirement leads to the conclusion that fault tolerance is not possible 

The second operational realization of quantum information involves specifying an ideal de- 
coding procedure. Such a decoding procedure adjoins I and (possibly) an ancilla system A to P, 
where I and A are in specified initial states 10)^ and |0)^. The total state space is determined by 
Hp © Til © H/K. The decoding operation is a unitary operator on this state space. After it is 
applied, the desired quantum information resides in I. The decoding -based view of quantum in- 
formation has been used successfully in analyses of fault-tolerant quantum architectures (see, for 
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example, i^). To connect the decoding-based realization of quantum information to subsystems, 
note that the decoding operation is, in effect, an isometry from an extended space TYp^ © TYj to 
H\ ® H\j. Here we have identified Hp^ with Hp ® |0)^ ® |0)^ and Hu with Ha ® Hp. Decoding- 
based realization is therefore equivalent to subsystem encoding in an extension of the physical 
state space, where the extension need not be physically meaningful. One can consider generaliz- 
ing decoding-based realizations by means of isometrics that provide the identification 

Hp^Hi®HT = Hv ^Hs® Hr, (2) 



which is obtained if the decoding operation also involves additional physical systems in unspeci- 
fied initial states. Although this is more general than subsystem encoding, most such isomorphisms 
do not result in quantum information that can be considered to be faithfully encoded in system P. 
We resolve this problem at the end of Section |n|by pairing decoding and encoding operations. 

A third approach to encoding of quantum information uses operators to characterize quantum 
information and makes it possible to give a reasonable definition of "faithful encoding of I in P". 
We give such a definition at the beginning of Sect. Inland prove that every such encoding is associ- 
ated with a subsystem. The intuition is that the states of I are characterized by the expectations of 
a linearly closed set of observables O. To ensure the correct dynamics of these states, the complex 
multiples of the observables should form an algebra f-isomorphic to B{H). Thus, one expects that 
faithful encoding of quantum information requires identifying a f-closed subalgebra A of B(Hp) 
that is isomorphic to B(H\). If this has been done, then the representation theory of f-closed al- 
gebras uniquely identifies a decomposition Hp = Hy ®Hs® Hr such that A consists of all the 
operators acting only on Hy. The results of Sect.Hllprovide support for the "subsystems principle" 
for realizing quantum information: 

The subsystems principle: Any faithful representation of quantum information in a physical 
system requires that at every point in time there are identifiable subsystems encoding the desired 
quantum information. 

The subsystem principle is powerful, but it is worth noting that it is sometimes convenient to 
use realizations of quantum information that do not satisfy this principle perfectly. For example, 
in optical quantum computing with "cat" states, it is convenient to represent the logical states of 
qubits by non-orthogonal coherent states Hldlllh . Another example is the study of "initialization- 
free" decoherence-free subsystems, where the probability amplitude of the encoded information 
may be less than 1 and the nature of the remaining amplitude must be taken into account 11211 . 



B. Protecting Quantum Information 



In most physical settings there are sources of errors that can affect encoded quantum informa- 
tion. Ideally, we would like exact knowledge of the error behavior of a physical system under all 
circumstances in which it is used. Since this knowledge is usually unavailable, one of a number 
of idealizing assumptions can be made. In the context of quantum channels, or when unwanted 
interactions are expected to have weak temporal correlations, we assume errors to be due to a 
known markovian process (in the continuous time setting) or a known quantum operation (in the 
discrete time setting). Both may be described by a collection of possible error events 8 = {Ei}. 
In general, the goal of quantum error control is to find quantum information subsystems for which 
the effects of the markovian process or quantum operation can be suppressed to the largest extent 
possible. Because the exact nature of the errors is usually not known, this goal is typically difficult 
to pursue. To make the task more tractable, we can consider only those errors that are expected to 
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be likely and look for subsystems that allow for "good" protection against such errors. We can then 
bound the effect of other errors by making estimates of their maximum probability (or amplitude) 
of occurrence. 

In this paper, we focus on subsystems that enable perfect protection against a fixed set of errors 
{Ei}i, with or without active intervention. Because of the linearity of quantum mechanics, perfect 
protection against the Ei implies perfect protection against any error in the linear span S of the Ei. 
Subsystems whose states are unaffected by the errors are known as "noiseless" or "decoherence 
free" subsystems [.13J and were introduced in |[l|,|7l] in the context of j-closed £ (or the j-closure 
of a non-f-closed £), in which case they can be characterized by irreducible representations of 
the commutant of £, which is the set of operators that commute with all members of £. In gen- 
eral, noiseless subsystems are not as easily characterized. In lfl4ll an explicit characterization of 
noiseless subsystems for any £ is obtained. This characterization is readily seen to be equivalent 
to the statement that the subsystem I' of the decomposition Hp = H-y ® Hs © is noiseless if 
and only if the restriction of the E^ to the subspace H\' ® acts as 1*^''^ ® E^^^\ Several equiv- 
alent characterizations for when £ is the span of the operation elements of a specific quantum 
operation were obtained in [0, IH]. These characterizations do not directly address the question 
of how one can computationally search for noiseless subsystems. A strategy for this search was 
offered in 0]. This strategy requires finding ^'-invariant subspaces and decomposing them into the 
canonical subsystems associated with the irreducible representations of a fixed-point algebra for a 
quantum operation whose operation elements span £. In Sect.|llllwe develop this strategy into an 
algorithm that does not require explicit constructions of algebras other than that generated by £. 
The mathematical structure of algebras over the complex numbers plays a crucial role. Interest- 
ingly, if there exists a quantum operation whose operational elements span £, then the algorithm 
simplifies substantially and is efficient in the dimension of the Hilbert space. Note that there is 
no a priori requirement that the likely errors included in £ be derived from a quantum operation. 
However, in most cases £ does satisfy this condition. To ensure that this condition holds, one can 
add I — A EjEi for a sufficiently small A, although the choice of spanning set Ei and A may 
affect the availability of large-dimensional noiseless subsystems. 

When no noiseless subsystem of sufficiently large dimension can be found, it is necessary to 
use active intervention to protect encoded quantum information. The idealized setting for active 
intervention involves alternating steps consisting of error events E^ and a quantum operation TZ that 
ensures that the errors do not affect the encoded information. An operation TZ with this property is 
known either as a "recovery" or as an "initialization" operation, depending on context. According 
to the subsystems principle, there must be two subsystems, one in which quantum information 
resides after error events but before TZ is applied, and another after TZ is applied. We call the first a 
"protectable" subsystem. The second is known as an "error-correcting" subsystem. Provided the 
encoded quantum information has been successfully protected, both subsystems are noiseless. The 
first is noiseless for the products EiRj, where the Rj are the operation elements of TZ, whereas the 
second is noiseless for the operators RjEi. Protectable subsystems are defined (but not named) 
in lU], where it was shown how to determine the protectable subsystem in the case where the 
error-correcting subsystem is a quantum code, that is, the cosubsystem is one-dimensional. Error- 
correcting subsystems are the main feature of operator quantum error correction 10, El] and directly 
generalize traditional error-correcting codes. 

Knowledge of the protectable subsystem and the error-correcting subsystem associated with a 
recovery/initialization operation and the relationship between the two helps us to understand how 
quantum information is stored at all times. An advantage of the protectable subsystem is that 
in many cases it is a simple extension of the error-correcting subsystem. That is, the former's 
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cosubsystem is a consistent extension of the latter's cosubsystem. As a result, the observables 
associated with the protectable subsystem induce the correct observables on the error-correcting 
subsystem. This implies that for the purpose of identifying the current value of the stored quantum 
information, it suffices to know the protectable subsystem, regardless of whether the last event 
was a recovery operation or an error. Examples of this situation are stabilizer codes with decoding 
algorithms based on syndrome extraction. It is readily verified that the associated protectable 
subsystem contains the stabilizer code as a subspace where the cosubsystem is in a particular state. 
In particular, this property holds for the stabilizer-based error-correcting subsystems identified 
in lllSh and used to simplify Shor's 9-qubit one-error-correcting code 1II6I1 . except that the error- 
correcting subsystem is defined by a subspace of the protectable subsystem's cosubsystem. In 
general, this relationship between protectable and error-correcting subsystems always holds if 
I G It becomes particularly useful in the context of fault-tolerant quantum computation, where 
the recovery operation and error events can no longer be easily separated. In this case the ideal 
error-correcting codes or subsystems associated with a scheme are typically not where quantum 
information resides. It resides in the protectable subsystems of the scheme. Note that in this 
setting it is usually the case that the subsystems containing quantum information vary in time. This 
happens, for example, when teleportation is used for error correction, when quantum information 
is stored in memory versus being actively manipulated, and in cluster-state-based schemes as part 
of the model fl?]. 

An advantage of error-correcting subsystems is that there are simple criteria and algorithms 
for determining whether there exists an associated recovery operation for which it becomes noise- 
less Not having to specify the recovery operation simplifies the search for subsystems 
suitable for protecting quantum information and makes it natural to talk about error-correcting 
subsystems without specifying the recovery operation. The same cannot be said for protectable 
subsystems. In Sect.|lV|we partially remedy this situation by reducing the problem of determining 
whether a subsystem is protectable to a number of other problems not involving the existence of a 
quantum operation. 

C. Conventions 

Capital letters in sans-serif font such as A, H, P are used to label quantum systems. The 
state space of a system A is determined by a Hilbert space, denoted by TYa- We label states 
according to the quantum system they belong to. For example, IV')^ is a pure state of A and p*^^^ is 
a density matrix for A. The tensor product symbol may be omitted in tensor products of labeled 
states and operators. We frequently consider instances of identical state spaces realized by and in 
different systems. We use primes to distinguish the different systems with identical state spaces. 
Thus, I, I' and I" are systems whose state spaces are identified via implicit isometrics, which are 
inner-product-reserving linear maps. In particular, a state of I is identified with the states li/j)^^ 

and [tp)^^^ of systems I' and I". One way to interpret this is to consider ^ as a symbol labeling 
a vector in an appropriate Hilbert space H and \'ip) l-ip)^ as the isometry identifying H and 
7ii. We use the equality symbol "=" not just to denote strict mathematical equality but also for 
identifying objects which are equal via an isomorphism. The isomorphisms involved are defined 
only implicitly, provided the meaning is clear. For a Hilbert space H, B{l-L) denotes the algebra of 
operators of H. U (H) denotes the group of unitary operators of H. In this work, all state spaces 
are finite dimensional. 
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II. FAITHFUL ENCODINGS OF QUANTUM INFORMATION 



To formalize the idea of "faithful encoding" we consider more general ways of encoding quan- 
tum information. A faithful encoding of I in P is a map D from density operators p on H\ to 
non-empty sets of density operators on Hp together with a map O from observables (hermitian 
operators) A of H\ to non-empty sets of observables of an extension Hq of Hp that satisfies three 
faithfulness requirements: 

1. Statics: For all a e D{p) and X E 0(A), 

tr{aX) =tr{pA). (3) 

This requirement ensures that we can identify the expectation values of faithfully encoded 
states. 

2. Unitary dynamics: For all a E D{p) and X E 0(^4), 

e-'^ae'^ E D(e-'^pe'^). (4) 

With this requirement satisfied, we can evolve the states using conventional quantum control 
so that the evolved states are consistent with the first requirement. 

For the next requirement, extend the domain of D to all positive semidefinite operators by defining 
D{0) = {0} and for p ^ 0, D{p) = tr(p)D(p/tr(p)). For an operator Z, let U{Z,X) be the 
projector onto the A-eigenspace of Z, or, equivalently, the projector onto the null space of Z — \. 
For A not in the spectrum of Z, the projector is 0. 

3. Measurement dynamics: For all a E D{p) and X E 0{A) and A real, 

n(X, A)(Tn(X, A) E D{I{{A, A)pn(v4, A)). (5) 

Faithful measurement dynamics ensures that projective measurements can be implemented 
correctly. 

The support of a positive semidefinite Hermitian operator p is the span of its non-zero eigen- 
value eigenvectors and is denoted by Supp(p). For a set of such operators D, Supp(Z}) is the span 
of the supports of the members of D. 

Theorem 1 If D and O are a faithful encoding of I in P, then one can identify a subsystem en- 
coding Hp = H\i ® Hs ® Hr such that for all p, D{p) has support in Hy ® Hs, and for all A, 
H\' ® Hs and Hr are invariant subspaces ofO{A), and 0{A) acts as A' on Hy ® 



The conclusion of the theorem does not hold if we assume only faithful statics and faithful 
unitary dynamics. For example, any irreducible representation of U{H\) leads to an encoding sat- 
isfying these two faithfulness requirements, and such representations that have dimension larger 
than N exist. For example, if I is a qubit, then any spin > 1/2 representation of SU{2) yields 
an encoding that lacks faithful measurement d yna mics. An other example is ensemble quantum 
computing with pure or pseudo-pure states 1 18, 1^ 23]. In the case of pseudo-pure states, faithful 
statics is only satisfied up to a scale. Nevertheless, quantum information is still associated with 
subsystems. It may be interesting to determine the nature of encodings satisfying only faithful 
statics (perhaps weakened to allow for an overall scale factor) and faithful unitary dynamics. Are 
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they always equivalent to a sum of subsystems transforming under distinct irreducible representa- 
tions of U (H\)l On the other hand, we conjecture that faithful statics and measurement dynamics 
imply faithful unitary dynamics. However, our proof of the theorem requires all three faithfulness 
properties. 

An apparently more general faithfulness property, "faithful interactions", requires that the en- 
coding of I in P behaves correctly in interactions with other idealized systems. Faithful interactions 
are needed if we use the encoded quantum information in a quantum information processing set- 
ting with multiple physical systems, each encoding quantum information in some way. Faithful 
measurement dynamics can be seen to be a special case of faithful interactions, and, according to 
the theorem, it implies faithful interactions in general. 

Proof of Thm.\I\ Let V be the linear sum of the supports of operators in D{p) for all p. Let V"*" 
be its orthogonal complement. By assumption, V C Hp. The proof proceeds in three stages. In 
the first, we show that the operators of 0{A) have V and V"*" as invariant subspaces. We can then 
redefine 0(^4) by restricting its operators to V. Hr is identified as V"*- fl Hp. We then show that 
0(^4) consists of exactly one operator and deduce that O extends to an algebra isomorphism when 
restricted to commuting subsets of observables. The underlying reason for this involves showing 
that the eigenspaces of 0{A) may be faithfully identified with eigenspaces of A. The first two 
stages of the proof do not require faithful unitary dynamics. The last stage involves analyzing 
SU (2) subgroups of U (H\) and corresponding subgroups of U (V) induced by O. Their action on 
eigenspaces of operators in the range of O implies the desired subsystem encoding. 

For an operator X, let Null(X) denote the null space of X. 

Lemma 2 Let X G 0{A). Then Null{X — A) fl V is non-empty if and only ifX is in the spectrum 
of A. Furthermore, V = {Null{X — A) fl V), and Null {X — A) fl V is the linear span of the 
supports of p G -D(cr) with Supp{a) C Null{A — A). 

Proof. Suppose that A is in the spectrum of A, and consider any p G D{Il{A, \)). By 
faithfulness of measurement dynamics, Il{X, \)pIl{X, \) G -D(n(y4, A)). By faithfulness of 
statics, n(X, A)pn(X, A) is not zero. Since the support of Il{X , \) pIl{X , \) is contained in 
Null(X — A) n V, this intersection is non-zero. Conversely, suppose that Null(X — A) fl V is 
non-empty. Then there exist a and p G D(a) such that the support of p is not orthogonal to 
Null(X - A). Thus n(X, \)pU{X, A) is not zero and is a member of D(Il{A, X)aU{A, A)). Be- 
cause D{0) = {0}, Il{A, A)crn(y4, A) is not zero. Hence Null(y4 — A) is non-zero, so that A is in 
the spectrum of A. 

To prove that V is spanned by the subspaces Null(X — A) fl V, we use the following sequence 
of inclusions: 



J2 (Null(X - A) n V) 



A 



EE E (Supp(n(x,A)pn(x,A))nv) 



A cr p:pGD((7) 



E E E(Supp(n(x,A)pn(x,A))nv) 



o" p:p&D{a) A 




(6) 



V, 



(7) 
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where in each expression, A ranges over the spectrum of A. The critical step in the sequence 
requires the inclusion 

J2 (Supp(n(x, \)pu{x, A)) n V) 3 Supp(p). (8) 

A 

To prove this inclusion, observe that Supp(n(X, X)pIl{X, A)) C V because n(X, X)pIl{X, A) e 
D{Il{A, \)pIl{A, A). Faithfulness of measurement dynamics implies that n(X, A)pn(X, A) = 
for A not in the spectrum of A. It then suffices to recall that for a complete set of orthogonal 
projectors P„ J2.Supp{P,pPi) = Supp(^,. P,pi^) D Supp(p). 

For the last claim of the lemma, let W be the set of density operators with support in 
Null(A -\).lf(reW, then U{A, \)aU{A, \) = a and for all A' ^ A, U{A, \')aU{A, A') = 0. 
Faithful measurement dynamics imply that for p G D(a), Supp(p) C Null(X — A). Thus, 
J2aew J2peD{a) Supp(p) C Null(X — A) n V. The following sequence of relationships proves 
the reverse inclusion: 

Null(X - A) n V = Null(X - A) n I E Yl Supp(p) 

c u{x, \)Iy H S"PP(/^) 

\ p:peD{a) 

= YY1 Supp(n(x,A)pn(x,A)) 

CT p:peD{a) 

^ E E supp(p) 

o- p:peD{U{A,X)crn{A,X)) 

= Y H ^"pp^/')' 

aeW p:p&D(a) 

where we have used the fact that for a projector IT and a positive semidefinite hermitian operator 

p, nsupp(p) = Supp(npn). ■ 



Corollary 3 Let X G 0{A). Then V and V"*" are invariant subspaces ofX. 

Proof. Because the eigenspaces of X are orthogonal. Lemma |2l implies that X can be block 
diagonalized with respect to an orthonormal basis whose first members span V. ■ 

LemmaEland Cor.|3limply that without loss of generality, we can assume that for all X G 0{A), 
X restricted to is 0. If not, replace every member of 0{A) with its restriction to V. This does 
not affect any of the faithfulness requirements. 

The last statement in Lemma [altogether with the assumption that X G 0(^1) has trivial action 
on implies that X's eigenspaces and eigenvalues are determined by A and the map D. It 
follows that 0{A) consists of exactly one operator. Thus, without loss of generality, we now take 
0{A) to be a function from observables of 1-L\ to observables of Tip. Lemma El also implies that 
inclusion relationships between eigenspaces of observables of 7Y| are preserved by O. 

Corollary 4 Suppose thatNull{A - Ai) C Null{B - A2). Then Null{0{A) - Ai) C Null{0{B) - 
A2). 
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If observables A and B commute, we can construct an observable C whose eigenspaces are 
the maximal common eigenspaces of A and B. By using the eigenspaces of 0{C) to derive the 
eigenspaces of 0{A) and of 0{B), we can see that 0{AB) = 0{A)0{B) and 0{aA + pB) = 
0{aA) + 0{PB), so that O preserves the algebraic structure of commuting sets of observables. 
Similarly, for any eigenbasis |Aj) of A, we can use an operator C with non-degenerate eigenvalues 
having the same eigenbasis to see that the spaces Supp(D(|Aj) (Aj|)) are a complete orthogonal 
decomposition of V into eigenspaces of A. From this it follows that 0(^4) is determined by the 
values of D on pure states. 

Forthelaststageof the proof of Thm.m we fix an orthonormal basis \i) ofH\. Lete^j = 

Xij = Cij + eji, Yij = —icij+ieji, C = Y^iiea and Vj = Null(0(C*)— i). Note that 2e~*'^^J''/^|i)| = 
g-jyij7r/2|^^ _ According to faithfulness of unitary dynamics, e~**^'^'^^-D(cr)e**^'^'^^ C 

D{e~'^"^ae^^). (For sets D and operators U, UD = {Ux:x E D}.) The inclusion is 
an equality because we also have Z}(e"*^cre'^) = e~**^(^)e*'^(^^L'(e"*^cre*^)e"*'^(^^e**^(^) C 
g-iO(A)^^^>)giO(A) rpj^-g jj^g g^^-gj. j-esults imply that e-'^^^^^^^'/'^Vi = Vj. By using Cor.gl 
with the eigenspaces of Xij, Yij and en + Cjj, we can see that the non-zero eigenspaces of 0{Xij) 
and 0{Yij) are contained in Vj © Vj. 

Because of the algebraic properties of O mentioned above and Xfj = en + Cjj, e~*'^(^'j)'^/2 = 
—iO{Xij). We can therefore fix an orthonormal basis \il) of Vj such that 0{XQi)\Ol) = \il) and 
0{XQi)\il) = \0l). The goal is to show that we can identify \il) with such that 0(^4) acts 

as the identity on the cosubsystem S. Note that the operators Xoj and Cjj generate the Lie algebra 
of U (Ti.\). Thus compositions of exponentials of the form e~*"^°J* or e~^^^^^ act transitively on the 
pure states of H\. It follows that for any \ip)^, Supp{D{\^jj)^\ilj\)) is an image of corresponding 

compositions of exponentials of the form e-«o(Xoi)t or e"*"^"" acting on Vq. Such compositions are 
completely determined by the basis Now 0(A) is determined by Supp (D(\ip)^\ijj\)), with \tp)^ 
ranging over eigenvectors of 0(A). Since for fixed /, O(Xoi) and O(ejj) act as they should on the 
states \kl), we have that 0(A) necessarily satisfies {kl\0(A)\kl) = \k\A\k)^, as desired. ■ 

We return to the issue of the relationship between decoding operations and faithful encodings. 
Decoding as defined in the introduction is the traditional way of identifying quantum information 
and generalizes recovery operations. A general form of the situation addressed by decoding in- 
volves one or more encoding isometrics Ci : H\ ^ Hp, one or more possible events Ei that are 
operators on Hp that may occur before we decode, and a decoding operation D that (after purifi- 
cation, if necessary) isometrically maps Hp into H\ ® Ha for some possibly composite system 
A. We say that {Q}, {Ej}, D preserve quantum information if for all DEjCi\tp)^ = \ip)^ 
for some unnormalized vector that does not depend on \ip)^. Here, which event Ei occurred 
is assumed to be unknown. We could consider the case where the decoding operation is chosen 
after the events and depends on partial knowledge of the events. However, by conditioning on the 
knowledge, we return to the situation just described. The only difference is that the subsystem 
associated with the situation may depend on the partial knowledge. To capture the case where 
quantum information is stored in error-correcting subsystems, let Ci be given by the isometrics 
identifying H\ with H\' ® 1^)^, where the range over any spanning set of Hs. 

In order to justify the subsystems principle, we prove the next theorem. 

Theorem 5 If {Ci}, {Ej}, D preserve quantum information, then there exists a subsystem encod- 
ing Hp = H\i ® Hs' © Hr such that for all i,j and EjCil^p)^ = ^ '^p 
= \'^\\(t^ij\' and \(j)ij)^ do not depend on 
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Proof. This follows from the fact that there exists a protectable subsystem associated with any 
quantum error-correcting code and associated recovery operation, which was proven in JH EUl . 
Alternatively, we could prove the theorem from Thm. [l]by defining D{\il))^ = {EjCi\tp)^} and 
0{A) by puUback of the appropriate operators via the decoding operator D. Here we give a direct 
proof. Let the \ be as required according to the definition of preserving quantum information. 
Let S be the set of vectors \(p)^ such that T-[\ ® |</))^ is contained in the range of D. Then S 
contains the Because the range of D is linearly closed, so is S. Define = S. Using the 

isometric properties of D, we can define a subsystem encoding Hp = H-v ® Ti-s' © Hr such that 
D(\'ip)^ 10) ,) = 1^) 10) • This subsystem encoding has the desired properties. ■ 



III. FINDING NOISELESS SUBSYSTEMS 

If ^ is a j-closed subalgebra of i?(7ip), the canonical decomposition of Hp is 

Hp = ^Hv,®Hs,®Hh, (10) 

i 

where operators A G ^ act as I^''^ ® S'.j(y4)''^'^ + 0*^'^'. For every operator of the form ^ . 1*^''^ ® 
Bi^^'^ + O^'^^ there exists an A G ^ with Si{A) = Bi. The H\. are noiseless subsystems for A. 
We also consider 7-^ r to be noiseless for A, but note that error operators in A have probability 
zero for states in this subspace. The tensor products and direct sums in the decomposition must 
be consistent with the Hilbert space's inner product. This is implicit in the construction and the 
identification via an isometry. 

Let £^ be a linearly closed set of error operators in B(Hp). For now, we do not assume that 
8 is the span of the operation elements of a quantum operation. Let Hy ® Hs © Hr define a 
subsystem encoding of I in P. Let 11 be the projector onto the support of I', Hy ® Hs C Hp. The 
subsystem is noiseless for S if and only if for all E E £, the restriction of E to Hy © acts 
as the identity on H\. Equivalently, for all E E S, EU = I*^'^ © S{Ey^\ It is straightforward 
to verify that if the subsystem is noiseless, then 11 projects onto an invariant subspace of S and 
£11 generates a f-closed subalgebra of operators acting on the support of 11 whose canonical 
decomposition contains noiseless subsystems with state space dimension at least N, the dimension 
of H\. Such noiseless subsystems are also noiseless for £. This leads to a strategy for finding 
noiseless subsystems with maximum dimensional Hy that is equivalent to the strategy proposed 
in 10]: 1. Pick an invariant subspace of £ and let 11 be its projector. 2. Determine the canonical 
decomposition of the f-closed algebra generated by US. The noiseless subsystems of this algebra 
are candidate noiseless subsystems for £. Our goal is to provide an explicit algorithm for finding 
suitable 11 and associated subsystems. The algorithm involves the decomposition of a matrix 
algebra, for which efficient algorithms are known, as we explain below. Note that in addition to 
the noiseless subsystems identified in this way, one can construct other noiseless subsystems as 
subsystems of already obtained noiseless subsystems, or by combining cosubsystems of identical 
dimensional noiseless subsystems with orthogonal supports. These constructions cannot yield 
larger dimensional noiseless subsystems, but they may generate ones with greater error tolerance 
or more efficiently controllable states. 

Let A be the algebra generated by £. Any noiseless subsystem for £^ is a noiseless subsystem 
for A. A is not necessarily f-closed. As a result, A does not have a canonical decomposition 
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of Tip as a direct sum of tensor products of Hilbert spaces. Nevertheless, we can identify a spe- 
cial subspace within which a similar decomposition is possible and where maximum dimensional 
noiseless subsystems may be found. This subspace is the span S of the irreducible subspaces of 
A. A subspace V of Tip is "irreducible" for A if it is ^-invariant, AV ^ and there is no non-zero 
.4-invariant proper subspace of V. Let Z be the null space of A. Both S and Z are invariant. 

Lemma 6 A maximum dimensional noiseless subsystem for A can be found in S or in Z. 

Note that Z is itself a noiseless subsystem. This subsystem is trivial in the sense that the 
probability of S-quoxs, is zero for any state in Z. This means that in a realistic setting, there 
must be operators acting on the system not included in S, and for Z to be at least approximately 
noiseless, they need to act as operators close to the identity when restricted to Z. 

Proof of Lemma^ Suppose that 7i|/ is a noiseless subsystem of Tip with cosubsystem Tis. 
Then V = Tiv ® Hs is invariant under A and for A e A, A acts as I^') (g) ^(A)^^^ on V. If for 
all A E A, S{AY^^ = 0, then V C Z and we are done. If not, then there exists a nontrivial 
irreducible subspace Si of Hs under the action of {S{A)A E A}. For each state \^)^,, l^jj)^^ ® Si is 
an irreducible representation for A. In particular, Si C S. Since H\' Si is also a noiseless 
subsystem, the proof is complete. ■ 

According to the theory of i?-modules, 5 is a module for A. The definition implies that it 
is semisimple, from which it follows that S = J2i '^i where the Si are irreducible and the sum 
is over independent subspaces, see [22], Chapter 9. Si and Sj are isomorphic with respect to 
the action of A if there exists an invertible linear map Uij from Si to Sj such that for \x) E Si, 
AUij\x) = UijA\x). The map Uij is said to "intertwine" Si and Sj. We can relabel the Si to 
form sets {Sik}i of isomorphic irreducible representations. For each k, let Vk be the span of 
the Sik and let [/q^^ be an intertwiner from Sok to Sjk- Choose a basis \iOk) of Sok and define 

\ijk) = uj^'^^liOk). Note that these vectors need not be orthogonal or normalized. Nevertheless, 
they define invertible linear maps from tensor products J7fc ® Sok to Vk via the linear extension of 
|j) ® \iOk) \ijk). The action of A E A with respect to this factorization is on Sq^ only. 

Lemma 7 A maximum dimensional noiseless subsystem for A in S can be found in one of the Vk- 

Proof. This follows from the argument given in the proof of Lemma|6l It suffices to observe 
that the irreducible representations ® Si are isomorphic for different ■ 

The main remaining problem in narrowing the search space for maximum dimensional noise- 
less subsystems is that the factorization of the Vk may fail to preserve the inner product. To 
simplify the notation, fix k and let V = Vk, Sq = Sok and J = Jk- Let U be an invertible linear 
map from J ® SqXoV that implements the above-mentioned factorization of V. Thus, for A G ^ 
and \x) E J ® 5o, AU\x) = ® -R(A))|x), where i? is a well-defined, irreducible represen- 
tation of A on iSq. Note that an irreducible representation of A on S^ is onto B{So) (Burnside's 
theorem). This implies that any noiseless subsystem of V must be associated with a subspace J'' 
of i7 such that the restriction of U io J' ® So has the property that there are linear operators W on 
J'' and V on Sq such that U{W ® V) is an isometry. Fortunately, in cases where A is generated 
by the operational elements of a quantum operation, we do not need to search for such subspaces. 

Lemma 8 If A is generated by the operational elements of a quantum operation, then there exist 
linear operators W on J and V on Sq such that U {W ® V) is unitary. 
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Proof. Let {Ei}i generate A, where the Ei are the operational elements of a quantum oper- 
ation O. By composing O with itself sufficiently many times, it is possible to obtain a quantum 
operation O' such that its operational elements span A. Thus, without loss of generality, assume 
that the span A and Y.,. E\Ei = I. We have Ei = U{I ® R{Ei))U'^. In order to continue, 
assume, without loss of generality, that J' ^ Sq = V. This can be done by means of any isometry 
between V and J' Sq. This implies that U is an invertible but not necessarily unitary linear map 
from J' Soto itself. We have 

U^^\l ® R{E,)^)U^U{1 ® R{E,))U^' = I, (11) 

i 

or, equivalently, 

^(I ® R{Eiy)U^U{l ® RiEi)) = U^U. (12) 

i 

This implies that for all positive semidefinite a on J', 

J2 R{Ei)hxj{{a ® l)U^U{a ® l))R{Ei) = tr j{{a ® l)WU{a ® I)), (13) 

i 

where trj is the partial trace over J'. Let R be the operation defined by R(X) = 
Y,i R{Ei)^XR{Ei). According to Eq.[ni for all positive semidefinite a, \Xj{{a ® l)WU (a ® I)) 
is a positive semidefinite fixed point of R. The spanning assumption on the Ei and irreducibility 
of So under R(A) imply that the R{Ei) span B(So). It follows that if p 7^ is positive semidef- 
inite and R(p) = p, then the support of p is So- It also implies that R has at most one positive 
fixed point (up to positive multiples): If p' is another one, then so is p — ep' for all e. Let e be the 
largest such that p — ep' is positive semidefinite. Then p — ep' is a fixed point with non-maximal 
support, which implies that it is 0. Let p be the unique trace 1 positive fixed point of R. Then, 
for all positive semidefinite a, tTj{{a ® T)WU{a ® I)) is a multiple of p. We can now deduce 
that WU = p' ® p for some strictly positive p'. Defining V = p~^/^ and W = p'~^/^ yields the 
lemma. ■ 

The above suggests the following strategy for finding maximum-dimensional noiseless subsys- 
tems: 1. Determine the span S of the irreducible subspaces of A. 2. Decompose S into a direct 
sum 0jXj of subspaces spanned by isomorphic irreducible subspaces. 3. For each Xj, let Ai be 
the restriction of A to and find the canonical decomposition for the j-closed algebra generated 
by Ai. This strategy will find maximum-dimensional noiseless subsystems provided that A is 
generated by the operational elements of a quantum operation. There are efficient algorithms for 
each step of this strategy; for a review, see 11231] . For completeness, we outline an algorithm that 
implements the strategy. 

To find S, consider the structure of A in more detail. If A does not contain I, replace A by 
A + CI. By doing so, the action of A on Z is no longer zero, but Z is still distinguishable 
from the other irreducible subspaces. Every one-dimensional subspace of Z is irreducible and not 
isomorphic to the irreducible subspaces of S. There exists a maximal chain of invariant subspaces 
= Vo C Vi C . . . C V„ = Hp such that the action of A induced on the quotients Vk+i/Vk is 
irreducible or zero. In a basis \ekj) of Hp where \ei^k+i)j) ^ V^+i \ (\ denotes set difference), the 
operators of A are block upper triangular. Let J' be the members of A that act as on each of these 
quotients. J' is known as the Jacobson radical of A. Let J\f be the null space of J', which is the 
set of vectors in the intersection of the null spaces of operators of J'. Then J\f is invariant (because 
JT" is a two-sided ideal) and S C Af (because S is invariant and the span of irreducible subspaces). 
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A fundamental property of J is that Aj J a semisimple algebra. Let Aj^ be the restriction of 
A to M . Then Aj^ is isomorphic to a quotient of A/ J, which implies that Af^ is semisimple. 
According to the representation theory of semisimple algebras, A/" is a semisimple ^^v^-module, 
which implies that N = S. Thus, to determine S, we can use an efficient algorithm for finding the 
Jacobson radical and then compute its null space. Decomposing M into independent irreducible 
subspaces can be done by means of an efficient algorithm for the decomposition of semisimple 
algebras over the complex numbers. A randomized algorithm can be based on the observation that 
Aj\f is isomorphic to a direct sum of complete matrix algebras Ak on the Squ, acting canonically on 
the (non-unitary) decomposition of S into a direct sum of products 1^ = Jk® <Sok- It follows that 
a random matrix in A// (with respect to a suitably chosen probability distribution) typically has 
generalized eigenspaces that generate (by multiplication by members of ^^v") exactly one of the 
invariant subspaces = Jk® <Sok- This yields the desired matrix algebras Ak- For each Ak, let 
Al be the f-closed algebra generated by Ak- The canonical factorization of Ik with respect to Al 
can also be obtained by a randomized algorithm. By construction, = Tii (g) 7-^2 (isometrically), 
with Al acting only on 7-^2- The eigenspaces of a randomly chosen Hermitian operator H2 in Al 
are typically of the form Hu \i) for an orthonormal basis of 'H2, where Hu = Hi, but with 
the isometry for making this identification not yet known. With high probability, these isometrics 
can be determined from an independently chosen second H2 by expressing H2 in an orthonormal 
basis whose i'th block of vectors is a basis of Hu Because H2 is a Kronecker product with 
identity action on Hi, the i,j block of H!^ must define an isometry between Hu and Hij (if it is 
nonzero). These isometrics must be consistent and induce the desired tensor product structure. 

Components of the algorithm of the previous paragraph not given explicitly include the gener- 
ation of an algebra from a set of matrices (this comes up in generating A from an error set and 
generating j-closed algebras from a given one) and various standard matrix manipulations such as 
matrix multiplication, eigenvalue and eigenspace determination, etc. We do not discuss the latter 
here. To generate the matrix algebra from a set of operators {E,i}, assume without loss of general- 
ity that the Ei are independent. Then iteratively, choose i, j and determine whether EiEj is in the 
span of the E^. If not, adjoin it to the set. Stop when for all i, j, EiEj is in the span of the E^. 

TV. PROTECTABLE SUBSYSTEMS 

As above, let {Ei}i be a set of error operators on Hp- Let Hp = Hy 0Hs(BHr be a subsystem 
encoding. The subsystem I' is "initialization protectable" (or "protectable" for short), if there 
exists a quantum operation with operation elements {Ri}i such that I' is noiseless for {EiRj}ij. 
The goal of this section is to reduce the problem of determining whether a given subsystem is 
protectable to the problem of searching for certain extremal error-detecting codes. We then reduce 
this problem to several linear algebra problems. 

Let be an orthonormal basis of Hs- For any state [tp) of H\' ® Hs C Hp, we define 

^(illV^) G H\' by the identity ® K)g = 1^)- V be the intersection of the inverse 

images of Hy ® Hs under the errors Ei. 

Lemma 9 With the definitions of the previous paragraphs, I' is protectable if and only if there 
exists a subspace V C V with the property that the maps Eij : \ip) ^j\Ei\ip) are proportional 
to a single isometry from V to Hy. 

Proof. For the "if" part of the lemma, we show that V is an error-correcting code for {Ei}i. 

We can reconstruct Ei on V from the Fij by the identity Eilip) = J^ii^ijli')) ® \j) ■ ^ be 

s 
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the isometry such that ^{j\Ei\i)) = aijU\ij). Then Ei\ij) = {U\ij)) ® Xl, «iib') • That V is an 
error-correcting code follows immediately. The operators Ri are given by U~^\i)^{i\. 

For the converse, we can use the subsystems principle (more specifically, Thm.|5l), according to 
which there must be a subsystem decomposition Hp = Hy (^Hj®Hq such that Ridip)^^ ® = 
l^)|„ ® I'Pij^ and Ei\'ip)^Jj)_^ = \'ip)^^\(pij)^. The desired subspace is given by ® |0)^ for any 
base state |0)^ of T. Note that the desired isometry is implicity defined via the two subsystem 
decompositions. ■ 

The maps Fij defined in the statement of Lemma|9lare well defined from V to'H\. Let M, N be 
the dimensions of V and H\, respectively. By choosing orthonormal bases of V and 

of H\i, the Fij are expressible as x M matrices (also denoted by Fij) with entries {Fij)ki = 
Without loss of generality, M > N, for otherwise the subsystem I' is clearly not 
protectable. The condition in Lemma|51can be seen to be equivalent to the requirement that there 
exists a unitary matrix V such that FijV contains a multiple of the N x N identity matrix as its first 
N X N block. The code V is spanned by the first N columns of V. This requirement is reminiscent 
of the familiar condition on the existence of an A^-dimensional error-detecting quantum code, 
according to which there must exist a unitary matrix W such that WEiW^ has a multiple of an 
N X N identity matrix as its first diagonal subblock. The protectability requirement can indeed 
be reduced to the existence of an error-detecting code. In particular, I' is protectable if and only 
if there exists an A^-dimensional error-detecting code for the operators {Fl-Fki}. Note that this is 
equivalent to requiring the existence of an A^-dimensional error-correcting code for the operators 
F^-, where Fl- is the square matrix obtained from Fij by expanding with rows of zeros. However, 

we do not have to consider all operators F-jFm. It suffices to find an A^-dimensional error-detecting 

code for operators of the form F-jFij and F-^F^ij), where tt is a cyclic permutation of the index 
pairs. 

We call subspaces V satisfying the condition in Lemma|51 "protecting" codes (for I'). There are 
several procedures that can be used to reduce the difficulty of the search for protecting codes. 

Lemma 10 All protecting codes are contained in the null space of the F in the linear span of 
{Fij}ij whose rank is strictly less than N. 

Proof. Let V be as specified in the paragraph before the statement of the lemma. If the rank 
of F is less than N, then the first N x N block of FV must be zero, from which the result follows. 

■ 

Let Gi, . . . , Gfc be X M matrices. We say that {Gi, . . . Gk] has maximal row rank if the 
span of the rows of the Gi has dimension kN . The next lemma generalizes Lemma flOl 

Lemma 11 Let Gi, . . . ,Gk be in the span of the Fij such that {Gi, . . . , Gk} does not have maxi- 
mal row rank, but for every k — 1 independent G[, . . . in the span of the Gi, {G[, . . . 
has maximal row rank. Then any protecting codes are contained in the intersection of the null 
spaces of the Gi. 

Proof. Let V be such that G,iV has an initial block proportional to the N x N identity matrix 
and V is spanned by the first columns of V. Suppose that V is not contained in the null space 
of some Gi. Then GiV's initial N x N block is not zero. The space Q of matrices G in the span 
of the Gj such that GV has an initial N x N zero block is A; — 1 -dimensional. Because the row 
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span of Gi is independent of the linear span TZ of the rows of the matrices in Q, the dimension of 
IZ is strictly less than [k — l)N, contradicting the assumption of the lemma. ■ 

LemmaE] means that in principle, the problem of finding V can be reduced to the case where 
each Fij has full rank and its row space is independent of the space spanned by the rows of the 
other Fki. In this case there are at most M /N independent matrices Fij. Unfortunately, we do not 
know of an efficient algorithm for checking the condition of Lemma[TT]that would enable reducing 
the problem to this case. Nevertheless, we can show that one can reduce to the case where there 
are at most M — 1 independent F^j. 

Lemma 12 For N > 1, if there are M or more independent Fij, then there exists a nonzero G in 
the span of the Fij such that G does not have full rank. 

Proof. Let {Gi]\^i be a basis of the linear span of the Fij. Let gl be the j'th row of Gj. If one 
of the gl is zero, we are done. Suppose / > M. Then the g] are dependent, so there is a non-trivial 
linear combination of the Gi with zero first row. Suppose / = M. Consider the matrices Aj whose 
i'th rows are the gf. Then there exists a non-zero linear combination aAi + jSA-z with determinant 
zero. Let x 7^ be in the null space of {aAi + 13A2Y . Then G = J2i ^i^i is not zero and the row 
vector y = (a, /3, . . .) satisfies yG = 0, so that G does not have full rank. ■ 

Note that the proof of the lemma contains an efficient algorithm for finding a non-full rank G. 
Let p^'^^^ be a density matrix on systems A and B. What states cr of A can be obtained by 
projecting B onto a pure state lip) ? The "pure cr -projection problem" for p^^^^ is to determine a 

state l^jj) such that ^{ip\p^^^^\t(j) = pa for some j> 7^ 0, if such a state exists. 

B B 

Theorem 13 If the span of the rows of the Fij is M -dimensional, the problem of determining 
whether subsystem I' is protectable is efficiently reducible to a pure 1-projection problem. 

If the rows of the Fij do not span the full space, then the protectability problem may be reduced 
to a generalization of the pure I-projection problem. However, in situations where the original 
error operators are associated with quantum operations, the Fjj 's do not have a common null space, 
even after the restrictions of the previous lemmas have been applied. Otherwise there would be 
states for which all Ei have zero probability. 

Proof. Let Gi,...,Gk be a basis for the linear space spanned by the Fij. We can 
choose an orthonormal basis of V such that in this basis, the matrices Gi have a block form 
[Gil, Gi2, . . . , Gii, 0, . . . , 0], where the Gij are x Nj matrices of full rank. We attempt to find the 
desired subspace V by choosing an orthonormal basis for V. Let X be the matrix whose columns 
are members of this basis. We wish to solve the k identities ail = GiX for X and a = {ai)i. We 
can write X in block form, X = [Xi ; . . . ; X^] , where Xi is A'^ x and the Xj are placed one 
above the other. The desired identities can be expanded as 

i 

ail = Y,G^,X,. (14) 

i=i 

The Xj can be eliminated by solving the equations in order. That is, from ail = GuXi we obtain 
«! = and Xi = if Xi 7^ N, and Xi = aiG^l otherwise. We write this as a linear constraint 
Li ■ a = and an identity Xi = aiGu, where Li may be "empty" (if Ni = N) and we set Gn 
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to be any left inverse of Gn. Once we have obtained Xj = amGjm and linear constraints 
Lja = for j < i, we can solve for Xi by substituing in Eg. (141 



-i-i j 

GiiXi = ail - ^ ^ amGijGjm- (15) 

j=l m=l 

The right hand side of this identity is a matrix Hi that depends linearly on a. GaXi = Hi can be 
solved if and only if the columns of Hi are in the column span of Gu. This condition yields a set 
of linear constraints Lja = 0. If the constraints are satisfied, then we can compute X^ = G'^Hi, 
where G'a is a left inverse of Gu. We can therefore define Gim by the identity X/ = J2m=i '^j^im- 
At the end of this process, the only free variables remaining are the aj, which must be chosen to 
satisfy the orthonormality constraint on X, J2i ^l^i = I- Expanding, we get 

^ ^ ajakGljGik = I, (16) 

i jk 

subject to Lia = for all i. If the linear constraints cannot be solved, we are done. Define p^^^^ 
by 

^(AB)^,^^^(A)t^(A)|^.^B^^|^ (IV) 

i jk 

where t is chosen so that tr(p('^^)) = 1. Any state \ip) in the subspace defined by Li\tjj) = (with 

B B 

Li defined with respect to the basis consisting of the |j) ) that solves the pure I-projection problem 

yields a solution for a by letting a, be a suitably scaled multiple of the coefficient of | j) in the 

solution. It follows that to complete the proof, it suffices to project p^^^^^ onto the subspace of 
B satisfying the linear constraints and renormalize the resulting positive semidefinite operator. 
This operator is a density matrix for which the pure I-projection problem is equivalent to the 
problem of whether I' is protectable. ■ 

The pure I-projection problem may be reduced to a problem of finding special matrices in a 
linear space of matrices. 

Theorem 14 The pure I-projection problem is polynomially equivalent to the problem of finding 
a matrix with orthonormal columns in a linear space of matrices. 

Proof. Consider the pure I-projection problem for p*^^^) . By purifing p^'^^^ with the addition of 
an environment E, we obtain a pure state 1^)^^^ whose reduced density matrix on AB is p^^^'> . The 
pure I-projection problem is now equivalent to the problem of finding \(f)) such that ^{(f)\\il)) is 

B ABE 

proportional to a maximally entangled state between A and E. Note that without loss of generality, 
the dimension of E is greater than that of A. Otherwise, the problem has no solution. We can 
expand everything in a basis for the different systems' Hilbert spaces: = ai\i)^, 1^)^^^ = 
X]iifc"^«i*:K) \j) 1^) • Let Mj be the matrix with coefficients {MAki = rriijk. The property that 

ABE 

^(^llV')^g^ is maximally entangled is equivalent to the property that J2j ^j^j has orthonormal 
columns. 

Given any set of matrices Mj we can reverse the reduction of the previous paragraph by setting 

Mj = tMj with t = 1/ {J2j tr(Mj^Mj)) to obtain a state such that its pure I-projection problem is 
equivalent to the problem of finding a matrix with orthonormal columns in the span of the . ■ 

Whether there is an efficient algorithm for finding a matrix with orthonormal columns in a 
linear space of matrices is an open question. 
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